Dynamic interference of photoelectrons produced by high-frequency laser pulses 
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The ionization of an atom by a high-frequency intense laser pulse, where the energy of a single- 
photon is sufficient to ionize the system, is investigated from first principles. It is shown that as 
a consequence of an AC Stark effect in the continuum, the energy of the photoelectron follows the 
envelope of the laser pulse. This is demonstrated to result in strong dynamic interference of the 
photoelectrons of the same kinetic energy emitted at different times. Numerically exact compu- 
tations on the hydrogen atom demonstrate that the dynamic interference spectacularly modifies 
the photoionization process and is prominently manifested in the photoelectron spectrum by the 
appearance of a distinct multi-peak pattern. The general theory is well approximated by explicit 
analytical expressions which allow for a transparent understanding of the discovered phenomena 
and for making predictions on the dependence of the measured spectrum on the pulse. 

PACS numbers: 33.20.Xx, 41.60.Cr, 82.50.Kx 



Dipolc transitions between bound electronic states of 
matter mediated by intense laser pulses are affected by a 
fundamental energy shift induced by the oscillating exter- 
nal field of the laser known as the dynamic or AC Stark 
effect 0-11. The AC Stark effect induced by optical laser 
fields is used in a myriad of contemporary experiments 
to hold and align molecules to shape potential en- 
ergy surfaces [lfj, to make rapid transient birefringence 
11 , IH , and to control quantum systems [H , 14 [ . What 
happens if the transition is to the continuum and not to 
a bound electronic state? Here, we investigate the pho- 
toionization of an atom by a high-frequency strong laser 
pulse, where the energy of a single-photon is sufficient 
to ionize the system. We demonstrate that the inter- 
play between the photoionization and the AC Stark shift 
results in a hitherto unrecognized dynamic interference 
of photoelectrons of the same kinetic energy emitted at 
different times. The effect is universal and prominently 
manifested in the photoelectron spectrum. 

In bound-bound electronic transitions accessible to op- 
tical lasers, the AC Stark effect detunes the laser fre- 
quency away from the field-free resonant frequency, an 
effect which has to be compensated in the presence of 
strong fields in order to have a transition. The situation 
is different in the case of ionization. The new genera- 
tion of light sources, like the free electron laser (FEL) at 
FLASH [lfj, and recentp rogress in high-order harmonic 
generation techniques (la . Il7| allow one to produce ultra- 
short laser pulses with single-photon energies well above 
the ionization threshold of any matter. This immediately 
raises the fundamental question of how the well-studied 
single-photon ionization process in weak fields will be 
modified in intense laser fields of high frequencies. We 
shall demonstrate here that single-photon ionization by 
high-frequency strong pulses is accompanied by dynamic 
interference caused by the AC Stark effect in the contin- 



In order to illustrate the new fundamental effect, we 
have chosen as an explicit example the photoionization 
of the hydrogen atom which is not affected by many- 
electron correlations and, thus, is amenable to exact com- 
putations. The AC Stark effect arises from the indirect 
coupling of nonresonant (nonessential) states that do not 
participate directly in the excitation process [l8[ . To ar- 
rive at an unequivocal description of the process, we have 
computed exactly the quantum motion of the whole dis- 
crete and continuum electron spectrum of hydrogen ex- 
posed to a strong pulse. For the sake of transparency of 
presentation we show below explicitly only the essential 
states (i.e., Is — > ep ionization) and employ the simpli- 
fying rotating wave approximation (RWA), which turns 
out to describe well the underlying physics. In the full 
calculations we have used also the nonessential states and 
discarded the RWA (see below and Appendix for details) . 

An atom initially in its ground electronic state |J) of 
energy of E] = (chosen as the origin of the energy scale) 
is ionized into the final electron continuum state \Fe) of 
energy IP + e (where IP = Ep — E] is the ionization 
potential and e is the kinetic energy of the photoelectron) 
by a linearly polarized laser pulse of carrier frequency 
to and pulse-shape function g(t). Following (lil- 24 1. the 
total wave function of the system as a function of time 
reads 



#(t) = ai(t)\I) + / dea e {t)\Fs)e 



(1) 



where ai(t) and a e (t) are the time-dependent amplitudes 
of the populations of the initial and final continuum 
states, the latter being dressed by the field (redefined 
by multiplying with the phase factor e lujt [2lj]). Insert- 
ing ^(t) into the time-dependent Schrodinger equation 
for the total Hamiltonian of the atom plus its interaction 
with the laser field, we obtain the following set of equa- 
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FIG. 1: Photoelectron spectra of the hydrogen atom exposed 
to Gaussian-shaped pulses of 30 fs (left panel) and 3 fs (right 
panel) durations with carrier frequency ui = 53.6057 eV com- 
puted exactly for different peak intensities. In weak- field 
single-photon ionization, the photoelectron line appears as a 
Gaussian peak centered around the energy of eo = w — IP = 
40 eV (indicated by the vertical broken line) and a width pro- 
vided by the pulse: FWHM of about 52 meV for the 30 fs 
pulse, and 0.52 eV for the 3 fs pulse. The presently computed 
strong pulse ionization spectra are, however, much broader 
(by several times) , shifted to higher energies from the central 
electron energy of eo = 40 eV, and posses distinct modula- 
tions of the electron intensity. It will be shown that these 
modulations are due to dynamic interference in the contin- 
uum. 



tions for the amplitudes (atomic units are used through- 
out) 



iai(t) 



ia e (t)= {KM g{t)ai(t) + (IP + e-uj)a B (t). (2b) 

£ o is the peak amplitude of the field £(t) = £q g(t) cos ut, 
and d e — (Fe\z\I) is the energy-dependent ionization 
dipole transition matrix element, which can be exactly 
computed for hydrogen. 

The calculations were performed for a Gaussian- 
shaped pulse, g(t) = e _t l T , and carrier frequency 
uj = 53.6057 eV, which is well above the ionization po- 
tential IP = 13.6057 eV. The photoelectron spectra com- 
puted for two pulse durations and different peak inten- 
sities (defined as / = £o/87ra) are depicted in Fig. Q] 
One can see that the effect of a strong laser pulse on 
the spectra is enormous. Instead of a single peak at 
£o = w — IP = 40 eV with the width of the pulse as 
naively expected from weak field calculations, the spectra 
are substantially shifted and spread, and posses distinct 
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FIG. 2: (Color online) Time evolution of the photoelectron 
spectrum computed exactly for the hydrogen atom exposed 
to a Gaussian-shaped pulse of 30 fs duration, of carrier fre- 
quency uj — 53.6057 eV and peak intensity of 5xl0 15 W/cm 2 . 
At very early times, the spectrum starts to develop symmet- 
rically around the central electron energy of eo = 40 eV. At 
later times, but before the maximum of the pulse arrives, the 
maximum of the spectrum moves to higher electron energies 
(see the spectrum computed at -10 fs). When the pulse ar- 
rives (the time fs corresponds to the pulse maximum), the 
spectrum is not any more symmetric, but still does not pos- 
sess intensity modulations. The first modulations start to 
show up in the spectrum only after the pulse maximum (see 
the spectrum computed at +10 fs) has arrived. 



modulations of the intensity. The effect is so big that it 
can easily be verified experimentally. 

To uncover the underlying physics, we show the time 
evolution of the photoelectron spectrum in Fig. [2] com- 
puted for the 30 fs pulse and peak intensity of 5 x 



de{\4£ Q } g{t)a e {t), (2a) 1Ql5 W { 



/ cm . At short times the spectrum develops sym- 
metrically around eo = 40 eV as expected for weak fields. 
Until the maximum of the pulse arrives the spectrum re- 
mains nearly symmetric, but its maximum continuously 
shifts to higher energies. Importantly, strong modula- 
tions of the intensity start to develop in the spectrum 
only after the maximum of the pulse has passed. This 
will be explained below. 

So far, we discussed the spectra computed exactly. 
Now, we would like to interpret the findings by introduc- 
ing approximations which allow us to arrive at explicit 
analytic expressions and proven to be surprisingly accu- 
rate. In the local approximation [25|, [26[, Eq. (l2"a|) reads 



idj(t)= (A-ir) g 2 (t)aj(t). 



(3) 



Here, the real term, Ag 2 (t), is nothing but the AC Stark 
shift of the energy of the ground state immersed into 
the dressed continuum. As usual, this shift follows the 
intensity envelope of the field 18j|, i.e. g 2 (t). Derivations 
of the expression for A can be found, e.g. in [l8j, and 
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in the present case it explicitly reads (without using the 
RWA) 



10 
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V \ de 



d £ £ 
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IP 



■ e — lj 



IP 



(4) 

where V stands for the principal value of the integral. 
In the RWA approximation the second term on the r.h.s. 
of Eq. (01 is neglected. The value of A is proportional 
to the peak intensity, A ~ £ 2 ~ /, and it is essentially 
nonzero owing to the energy dependence of the dipole 
matrix element d e . 

The imaginary term — ^Tg 2 (t) in Eq. (|3|) describes the 
losses of the population of the ground state by the ion- 
ization into all final electron continuum states \Fe). T 
can be obtained from A in Eq. flU). The explicit expres- 
sion for r was obtained in 21] in the local and RWA 
approximations and in the present case it reads 



2tt IKo^oI 2 



(5) 



Eq. ([5]) is also valid without using the RWA, since the 
second term in the energy integral (|4]) possesses no pole. 
The decay rate T is also responsible for the additional 
broadening of the photoclcctron peaks beyond the width 
given by the pulse intensity envelope. 
Equation © can be solved explicitly 



ai(t) 



=( _ iA -r/2)j(t) 



(6) 



where J{t) = j t _ oo g 2 {t')dt' grows smoothly with time. 
We find that the result (JSJ) is in excellent agreement with 
the exact calculations. It also simplifies greatly the com- 
putation of the spectrum. As it stands, a F (t) in Eq. (|2b[) 
can be expressed as an integral of a/ (t) [IjJ, |2l| 



a e (t) = -i (KM e~ l5t [ 

J — i 



g(t')ai{t')e ldt dt' , (7) 



where we introduced the abbreviation S = IP + e — lu = 
e — £o, which is the electron energy detuning from £o. 
Using now the explicit expression (|6]) makes the compu- 
tations of a e {i) and of the spectrum a(e) = |a £ (+oo)| 2 
rather straightforward. 

The photoelectron spectrum computed in the local ap- 
proximation is depicted in the upper panel of Fig. [3] 
The AC Stark shift A = 0.26 eV and the decay rate 
T = 0.044 eV computed via the above equations (j4j and 
([5]) have been utilized in the numerical calculations. It is 
seen that the spectrum calculated in the local approxima- 
tion (solid curve) describes very well the exact spectrum 
(open circles). 

Due to the success of the local approximation, we are 
able to gain deeper insight into the origin of the inter- 
ference effects leading to the intensity modulation in the 
spectrum by evaluating analytically a £ {t). For this pur- 
pose we notice that the integrand in Eq. ([7]) contains 
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FIG. 3: (Color online) Upper panel: The photoelectron spec- 
trum of the hydrogen atom exposed to a Gaussian-shaped 
pulse of 30 fs duration, carrier frequency uj = 53.6057 eV 
and peak intensity 5 x 10 15 W/cm 2 computed exactly (open 
circles) compared with the spectrum computed in the local 
approximation (solid curve) with the quantities F = 0.044 eV 
and A = 0.26 eV. Lower panel: The spectrum obtained 
in the stationary phase approximation via Eq. (JSJ) with the 
same parameters (solid curve). The two individual con- 
tributions to the spectrum f rom the negative and positive 
times t s = ±ti(e) = ±r i/ln[A/(£ — £o)]/2 describing in 
Eq. (|8} the contributions when the pulse arrives and ex- 
pires, respectively, are shown by broken curves. The en- 
ergy dependent interference term in Eq. ((8} is proportional 
to cos {-2(e - e )ti - A[J(-ti) - J{+ti)] - 7r/2} . In the en- 
ergy interval e— eo £ [0, A] for which a stationary phase can be 
found, its argument accumulates a phase of AJ(+oo), which 
for a Gaussian pulse is equal to At^tt/2. For the present 
parameters, the interference term accumulates a phase of 
14.85 ~ 2.36 -27T rad, and, thus, the spectrum exhibits around 
two and a half oscillations of the intensity profile as seen in 
the figure. 



a rapidly oscillating factor exp [iSt' — iA J{t')\ which is 
multiplied by a smoothly varying function. The main 
contributions to the integral stem from the times at 
which the phase $(i) — St ~ AJ(t) is stationary [27| . 
i.e. 4> = 0. The resulting stationary phase condition, 
Ag 2 (t s ) = S, has a transparent physical meaning. It de- 
fines the time t s (e) at which the energy of the ground 
state, continuously shifted by the AC Stark effect, moves 
across the energy position e of the dressed continuum 
state (S = e — eq). For any pulse there are at least 
two stationary points for each value of e: one, ii(e), 
when the pulse is growing, and another, £2(^)7 when it 
decreases. For a Gaussian pulse there are exactly two 
times, ti(e) = — t 2 (e). 

By collecting in the integral ([7]) these two main contri- 
butions at t s — ±ii(e), we obtain the following explicit 
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FIG. 4: (Color online) The intense laser pulse of sufficient 
high frequency to ionize the system by one-photon absorption 
induces a time-dependent AC Stark shift of the ground state 
of the system. Because of this dynamic shift, the photoelec- 
tron emitted along the pulse envelope has at every moment t 
predominantly a specific kinetic energy e which is determined 
by that t (at a given time t, the phase contributing to the 
integrand in Eq. (0 is stationary at e = so + <? 2 (t)A, for 
details see text). Owing to the fact that the pulse envelope 
first grows and then falls, there are for a Gaussian pulse ex- 
actly two times (indicated as t~ and t + in the figure) at which 
the emitted electron wave has the same energy e. These two 
waves emitted with a time delay with respect to each other 
interfere. This dynamic interference gives rise to the strongly 
modulated intensity distribution of the photoelectron spec- 
trum seen in Figs. [TJ3] The dynamic interference can be con- 
trolled by the intensity, duration, carrier frequency and shape 
of the pulse. A pulse envelope with more than one maximum 
leads to more than two waves which interfere and higher laser 
intensity to larger AC Stark shift A and to more oscillations 
of the spectrum. 



approximate expression for the spectrum 
d £ £o 



a(e) 



r/2j(t„) e i[Hts)±^ 



t„=±ti(e) 



The additional phase factors ± ? result from higher terms 
in the expansion of the phase $(i) around the stationary 
points ±ti(e) computed for the Gaussian pulse. The pho- 
toelectron spectrum evaluated via Eq. © is depicted in 
the lower panel of Fig. [3] (solid curve) . It is illuminat- 
ing to see that an explicit simple expression reproduces 
nicely the exact spectrum (open circles in the upper panel 
of the figure) which is complicated to compute. The indi- 
vidual contributions to the spectrum of the two terms in 
the expression ([5]) are smooth and rather boring (broken 
curves in the lower panel of the figure). 

Eq. allows one to uncover the physical origin of 
the strong modulations in the photoelectron spectrum. 
These are the results of the coherent superposition of 
two photoelectron waves emitted with the same kinetic 
energy at two different times. A cartoon visualizing this 
dynamic interference is shown in Fig. 21 This new physi- 



cal effect causes enormous qualitative changes in the pho- 
toelectron spectrum which can easily be verified experi- 
mentally. Importantly, the interferences vanish if the AC 
Stark shift in the continuum vanishes (i.e., if A vanishes). 
Then, there are no stationary points along the pulse. We 
would like to stress that the dynamical interference can 
be controlled by the shape of the pulse envelope and is 
influenced by the choice of the carrier frequency. We are 
convinced that dynamic interference is an universal ef- 
fect in ionization processes by high-frequency strong laser 
pulses and its experimental investigation is feasible. In 
the following we address these issues in some detail. 

The predicted phenomenon also persists in other sys- 
tems with more electrons. In the hydrogen atom dis- 
cussed above the residual ion does not posses electrons 
and, therefore, the dynamic interference is governed by 
the AC Stark shift g 2 (t)A which is a property of the 
ground state in the field. If the ion produced by pho- 
toionization does posses electrons, then this ion is obvi- 
ously also subject to an AC Stark shift g 2 {t)A[ on . In this 
case, the dynamic interference pattern in the spectrum 
is governed by the total AC Stark shift g 2 (t)ATot where 
Axot = A — Aj on [28[. As the ion is usually harder 
to ionize than the neutral system, one can expect that 
Axot will be substantial even if both A and A/ ora have 
the same sign. Moreover, depending on the photon en- 
ergy ui, the shifts A and Aj on can have different signs, 
enhancing interference effects (sec an explicit example in 
[2^). In many-electron systems, several electronic shells 
can be subject to the ionization, and different ionic states 
can be produced. Each ionic continuum will contribute 
to the shift of the ground electronic state leading all to- 
gether to a final unique value of A at the carrier frequency 
u>. In general, different ionic states \F a ) will experience 
different AC Stark shifts g 2 (t)Af on . The dynamic inter- 
ference pattern in each of the partial electron spectra is 
rg\ determined by different total shifts A£ ot = A — A" ora 
pjj . and one can expect dissimilar dynamic interference 
patterns in different partial photoelectron spectra. Ex- 
perimental investigations of dynamic interference will be 
possible, in particular if the studied ionic thresholds are 
not too close to each other, i.e., if the resulting partial 
spectra do not overlap. The AC Stark shift scales also 
with the probability to ionize the system. As photoion- 
ization cross sections usually decrease with the photon 
energy, the shift will be larger for frequencies uj which 
are not too far from the ionization threshold. In general, 
one can expect particularly large shifts, and hence more 
pronounced interference patterns, in systems with cross 
sections which vary strongly with the photon energy, e.g., 
in the vicinity of a shape resonance. 

To demonstrate the feasibility of experiments on dy- 
namic interference, we discuss the present status of FEL 



facilities. The currently operating FEL at FLASH 15] 
generates soft x-ray radiation in the energy range of 26 
180 eV with flux of about 10 13 photons per pulse with 
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durations of 10-70 fs [29]. Combined with appropri- 
ate focusing optics, peak irradiance levels of more than 
10 16 W/cm 2 can be achieved at present (H3|. The SASE 
FEL sources do not produce so far monochromatic radi- 
ation and their pulses consist of many spikes with ran- 
domly fluctuating properties. The impact of these prob- 
lems on the experimental output has been studied, e.g., 
inRef. [Hj]. Nevertheless, the pulse properties at FLASH 
are continuously improving, and an averagetemporal co- 
herence of 3-6 fs can already be achieved [32j. Moreover, 
several new techniques have been proposed for improve- 
such as the selection of a single ultrashort 
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ment 

radiation spike [Hf . In addition, the recently tested FEL 
facility FERMI at Elettra [35| operates in single-mode, 
providing coherent pulses (without spikes) with durations 
of 30 to 100 fs in the photon energy range of 12 to 413 
eV, and it is expected to provide an unprecedented flux of 
about 10 14 photons per pulse. Concluding, the pulse du- 
rations, photon energies, temporal coherence, and peak 
intensities necessary for experiments on dynamic inter- 
ference are available at present. 

In conclusion, intense high-frequency pulses give rise 
to a strong AC Stark effect in the continuum, which, in 
turn, modifies the ionization process dramatically. The 
effects are universal and can be verified experimentally 
by current techniques. Whenever matter is exposed to 
intense pulses, like in FEL experiments, ionization takes 
place, and the AC Stark effects in the continuum and 
the resulting dynamic interference will come into play. 
They will be essential for understanding the experimental 
results. The effects are thus both fundamental for light- 
matter interaction and of practical importance for the 
analysis of FEL experiments. 

We thank A.I. Kuleff for many valuable discussions. 

Appendix.-In this appendix we discuss how the full cal- 
culations were done. We have numerically solved exactly 
the system of equations which accounts for the interac- 
tion of all essential and nonessential electronic states of 
hydrogen with the laser pulse without employing the ro- 
tating wave approximation. For the peak intensities con- 
sidered here, the effect essentially stems from the quan- 
tum motion of the electron eigenstates restricted to ns 
and nd Rydbergs with n < 6, mp Rydbergs with m < 16, 
and, of course, ep continuum with e < 100 eV. Let us 
designate the time-dependent amplitudes for the popu- 
lation of the even ns/nti-states by a e (t), and for the 
odd mp/ep— states by a (t). The system of equations of 
motion of these states in the presence of the pulse reads 



; e (t) = £ e a e (*)+E (f ^) 9(f) + Oo(*) 

odd ^ ' 

i„(t) = E a (t) + ]T (^A g(t) (e*** + e"*") a e (t) 

piifn. ^ ' 



spectrum was discretized as \Fe) = \Fsk) ■ V 'Ae*, where 
\Fsk) represents the state in the continuum interval of 
Aek centered around the electron energy Ek ■ The energy 
steps Aefc were chosen linear with respect to y/e, and 
concentrated symmetrically around the central electron 
energy e . The convergence of the solution with respect 
to the electron eigenstates included, of the integration 
over time, and of the discrete representation of the 
continuous energy have been ensured. 



with boundary conditions ai s (— oo) = 1 and 
ai^is(^oo) = for all other amplitudes. The continuous 
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SEVERAL IONIZATION STEPS AND THRESHOLDS 



In the main manuscript we report on the theory of photoelectron spectra produced by laser 
pulses and their analysis. Numerically exact results were presented for the hydrogen atom 
which demonstrate the effect of dynamical interference beyond doubt. For transparency 
of presentation the equations were kept as needed for the hydrogen atom. Since this atom 
possesses only one electron, the produced ion cannot be further ionized. In this supplemental 
material we show how the present theory can be applied to systems with more electrons. 

We start our discussion by first reminding the ansatz for the total wave function in the 
manuscript (Eq. (1)) for the hydrogen atom 



where aj(t) and a e (t) are the time-dependent amplitudes of the populations of the initial 
and final continuum states. When plugged into the Schrodinger equation, these populations 
could be expressed as (see Eqs. (2b) and (3) in the manuscript and discussion therein) 



As explicitly shown in the manuscript, the quantum motion of the nonessential electronic 
states results in the time-dependent AC Stark energy shift, Ag 2 (t), of the ground electronic 
state (see Eq. (IS2al) ). Because of this dynamic shift, two photoelectron waves emitted at 
rising and falling sides of the pulse superimpose, resulting in the dynamic interference effect. 
The interference pattern in the photoelectron spectrum depends substantially on the shift 
A: The larger the AC Stark shift the stronger is the effect and the more oscillations appear 
in the spectrum (see the caption of Fig. 3 in the manuscript). 

In the presence of a strong pulse, a system with more electrons can undergo further 
ionization. In order to incorporate this fact in the theory, one has to include the doubly- 
ionized electron continuum states in the ansatz (jSTl) for the total wave function. This is 
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iai{t) = (A-fr) g 2 (t)aj(t), 



(S2a) 



ia £ {t) = {\d £ £ } g(t) aj(t) + {IP + e - u) a £ (t). 
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straightforwardly done. Following Refs. [S1-S6], the total wave function of the system as a 
function of time now reads 

ty(t) = ai {t)\I) + J dea £ (t)\Fe)e- iuJt + J J dede' b ££ ,{t)\Gee')e- 2iu)t . (S3) 

Here, b EE r(t) is the time-dependent amplitude for the population of the doubly-ionized elec- 
tron continuum state \Gee') which is already dressed by the field (i.e., by the energy of two 
absorbed photons 2u), and e' is the kinetic energy of the secondary electron emitted by 
the ionization of the ion. Inserting \l/(t) into the time-dependent Schrodinger equation for 
the total Hamiltonian of the atom plus its interaction with the laser field, and implying the 
rotating wave [S7] and local [S8,S9] approximations, we obtain the following set of equations 
for the amplitudes (atomic units are used throughout) 

taj(t)= (A - |r) 9 2 (t) a/(t), (S4a) 
ia £ (t) = {|4£ } g{t) aj{t) + (iP + [ Ajon - |r /on ] Q 2 (t) + s - w) a e (t), (S4b) 
ike>(t) = {\d e ,Ev} 9® a e(t) + (DIP + e + e' - 2w) 6 ee ,(t). (S4c) 
The last equation, Eq. ( 1S4cj) . describes the population of the doubly- ionized continuum 



states \Gee') by the subsequent ionization of the ion as a function of time. In this equation, 
the quantity DIP = Eq — Ej is the double-ionization potential of the system, and d e > is the 
energy- dependent ionization dipole transition matrix element of the ion. 

By comparing Eq. (1S4bl) with Eq. ( 1S2bl) . one can see that the ionization potential is 
subject to a time-dependent AC Stark shift A Ion g 2 (t). More precisely, the ionic threshold 
is shifted by the AC Stark effect. If the energy of a single-photon is large enough to ionize 
the ion, this shift is induced by the quantum motion of the nonessential doubly-ionized 
continuum states. If the photon energy u is insufficiently large to further ionize the ion, this 
shift also exists and is induced by the indirect coupling to the discrete and also continuum 
spectrum of the ionized system. In the case where the further ionization of the ion is possible, 
the energy of the ionic state is also augmented by the imaginary term — ^T Ion g 2 (t) seen in 
Eq. (]S4bj) . This term is responsible for the losses of the population of the ion by further 



ionization into all final doubly-ionized continuum states \Gee'). The situation is similar 
to that discussed in Eq. (jS2bJ) (or equivalently Eq. (3) in the main manuscript) where the 
ground state population decreases by ionisation losses. The expressions for A Ion and T Ion 
are analogous to those of the neutral system given in the manuscript, but for the ion. 
Let us now discuss the consequences of the above and compare Eqs. flS2a|S2bj) and 



Eqs. (IS4a|IS4bp . Because of the AC Stark shift A Ion g 2 {t) of the ionic state, the dynamic 



interference pattern in the photoelectron spectrum is determined now by the total AC Stark 
shift AT 0t = A — Aj on . As the ion is usually harder to ionize than the neutral system, one 
can expect that Axot will be substantial even if both A and A/ on have the same sign. More- 
over, depending on the photon energy, the shifts A and A/ on can have different signs. This 
may enhance the dynamic interference effect in the photoelectron spectrum. This situation 
can indeed be achieved systematically by choosing carrier frequencies below the ionization 
potential of the ion. In this case, a nonresonant coupling of the ionization threshold with 
ionic Rydberg states may result in a negative AC Stark shift Aj on of the ion. 

Until now we have discussed the situation of a single ionization threshold. A many- 
electron system generally exhibits several open ionization thresholds, and different ionic 
states can be produced by the photoionization. It is straightforward to incorporate this fact 
in the theory. One has just to sum over the final ionic states. The ansatz for the total wave 
function ( IS3l) now reads 

tt(t) = a/(t)|i}+5] J de a a £a {t)\F a e a )e- iujt + Y, J J de a de^b £a£ ^t)\G^ a e^e- 2 ^ , (S5) 

a a/3 

where the sums run over several final ionic continuum states \F a e a ) and also doubly-ionized 
continuum states \GpE a e'p). The derivation of the equation of motion is in analogy to those 
for a single ionization threshold in Eqs. (1S4I) above. Obviously, each of the ionic continuum 
states \F a e a ) will contribute by a partial shift A a to the unique AC Stark energy shift A of 
the ground electronic state: 

A = ^A a . (S6) 

a 

In addition, each ionic state will experience an AC Stark shift Aj onl which is a sum over the 



partial shifts A^ n associated with a given doubly-ionized continuum state \Gpe a e'a): 



In general, different ionic states will experience different AC Stark shifts A" on . The 
dynamic interference pattern in each of the partial photoelectron spectra with photoelectrons 
of kinetic energy e a will be determined by a different total shift Aj. of = A — Af on , and 
one can expect dissimilar interference pattern in different partial photoelectron spectra. 
Experimental investigations of systems with several open ionization thresholds would provide 
important information on the AC Stark shifts of different ionic states. The investigations 
will, of course, be possible if the ionic thresholds are not too close to each other such that 
the studied partial photoelectron spectra do not overlap. 

In conclusion, in many-electron systems each partial photoelectron spectrum a will ex- 
hibit its own dynamic interference pattern dictated by A^ ot = A — Af on similar to the 
hydrogen atom case studied explicitly in the main manuscript. 

DYNAMIC INTERFERENCE IN He 

We have exemplified our argumentation by as accurate as possible numerical calculations 
on the He atom. The calculation of the dipole transition matrix elements in the neutral and 
ionized He were performed in the Hartree-Fock approximation (which is, of course, exact for 
the He + ion). The photoelectron spectra of He computed at two different intensities of the 
pulse with 30 fs duration are depicted in Fig. [T] of this supplemental material. The photon 
energy u = 39.587 eV is above the ionization potential IP = 24.587 eV of the neutral He 
atom, but below the ionization potential IPi on = 54.418 eV of the He + ion, and the ion 
cannot be ionized. This frequency is also just below the first 2s/2p-Rydberg states of the 
ion (40.813 eV). This results in a negative AC Stark shift Aj on < for He + (ls 1 ) state (ionic 
threshold), which stems mainly from the indirect coupling with the He + (2p 1 ) state via the 
field. The shift A of the neutral He is positive and hence the resulting total shift in the 
photoelectron spectrum A^ot — A — Aj on is larger than that for neutral He. For the intensity 
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FIG. 1: (Color online) Photoelectron spectra of He exposed to a Gaussian-shaped pulse of 30 fs 
duration with carrier frequency u) = 39.587 eV computed for two peak intensities (indicated in the 
legend). The photon energy u is above the ionization potential of the neutral atom, but below the 
ionization potential of the He + ion, and the ion cannot be ionized. This frequency is just below the 
first 2s/2p-Rydberg states of the ion (40.813 eV), leading to a negative AC Stark shift Aj on < 
for the ion He + (ls 1 ) state. The shift A of the neutral He is positive, and, hence, the resulting total 
shift in the photoelectron spectrum Axot = A — A/ on (indicated in the legend for each intensity) 
is positive and enhanced. 

of 2.5 x 10 14 W/cm 2 the total shift amounts to A Tot = 0.02 - (-0.14) = 0.16 eV and for 
the intensity of 5 x 10 14 W/cm 2 to A Tot = 0.04 - (-0.26) = 0.30 eV. As a consequence, 
the computed electron spectra are shifted to higher energy from the photoelectron energy of 
6q = (jO — IP = 15 eV (indicated by vertical dotted line) one would have at weak intensity. 
One can see from this figure that the computed photoelectron spectra of the neutral He 
exhibit distinct multipeak patterns, which are due to the dynamic interference effect. We 
also notice that in He relatively low peak intensities of a few 10 14 W/cm 2 are sufficient to 
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arrive at substantial effects which can be verified experimentally by current FEL sources. 
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